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Abstract 

We propose a concrete model which exhibits ferromagnetism and electron-pair con- 
densation simultaneously. The model is defined on two chains and consists of the 
electron hopping term, the on-site Coulomb repulsion, and a ferromagnetic interaction 
which describes ferromagnetic coupling between two electrons, one on a bond in a chain 
and the other on a site in the other chain. It is rigorously shown that the model has 
fullypolarized ferromagnetic pairing ground states. The higher dimensional version of 
the model is also presented. 
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1 Introduction 



Recently, UGe2 PQ, URhGe [2], and UCoGe [3] were discovered to exhibit ferromagnetic 
superconductivity. Experimental results suggest that the same electrons are contributing 
to both ferromagnetism and superconductivity, and thus, in the superconducting phase of 
these materials, the electrons are considered to condensate into a spin-triplet pair state 
unlike usual low temperature superconductors, in which electrons are forming non-magnetic 
spin-singlet pairs. Microscopic explanation of this phenomenon is a challenge in condensed 
matter physics, but the problem is rather subtle and difficult since we have to treat spin- 
rotation symmetry breaking and electron-pair condensation simultaneously. In fact, the 
mechanism for ferromagnetism alone in itinerant electron systems has not yet been fully 
understood, although there have been some rigorous developments 01 El El El El E] in the 
Hubbard model, which is one of the simplest models of itinerant electron systems. As for 
electron-pair condensation, some models have been shown to have pairing ground states [TQl 
ITT1 [T2| [T3l [HI H5j, but little is known about spin-triplet pairing ground states. 

The Coulomb repulsion between electrons combined with the Pauli exclusion principle 
can generate ferromagnetism in itinerant electron systems [U El El [71 El [9] , but it is unclear 
at the present time that the same can also lead to ferromagnetic superconductivity. Since 
the electrons form spin-triplet pairs, it is expected that there appears a kind of effective 
ferromagnetic interaction between electrons in the ferromagnetic superconductors. Thus 
a model containing ferromagnetic interactions as well as the Coulomb repulsion between 
electrons will be suitable for an investigation as a first step toward understanding microscopic 
mechanisms for ferromagnetic superconductivity. 

In this paper we propose a model which has fullypolarized ferromagnetic pairing ground 
states. The model consists of the electron hopping term, the on-site Coulomb repulsion and 
a short-range ferromagnetic interaction term. It is possible to consider the model in any 
dimension, but for simplicity we mainly concentrate on a one-dimensional version here. In 
the one-dimensional case, our model is defined on two chains. The electrons can hop along 
the chain direction, feeling on-site repulsion, and furthermore the two electrons, one on a 
bond in a chain and the other on a site in the other chain, feel a ferromagnetic interaction (see 
Fig. [I]). It is shown in our model that, owing to the ferromagnetic interaction introduced here, 
the electrons form spin-triplet pairs, and then the fullypolarized ferromagnetic pairing state 
is selected by the on-site Coulomb repulsion as the unique ground state (up to degeneracy 
due to spin- rotation symmetry). It is also shown that the ground state exhibits off-diagonal 
long-range order associated with local spin-triplet pairs. 

This paper is organized as follows. In the next section we set up the model and state the 
main results. In section El we prove the main results. In section H] we make some remarks 
on our model and ferromagnetic pairing states. We also present models with anisotropic 
spin-interactions and show that the models exhibit other kinds of spin-1 electron-pair con- 
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Figure 1: The image of the ferromagnetic interaction considered here. The solid circles 
represent lattice sites, and thick arrows represent electron spins. The electron on a bond in 
a chain favors to couple ferromagnetically with the one on a site in the other chain. 

densation. In the final section we briefly comment on the higher dimensional version of the 
model. 



2 Definition and Main Results 

We start by defining a lattice of our model. Let Ai and A2 be linear chains with L sites. We 
label sites in Ai and A 2 by integers, and half integers, respectively, as Ai = {0, 1, . . . , L — 1} 
and A 2 = {1/2, 3/2, . . . , L — 1/2}. It is assumed that L is a positive odd integer, and periodic 
boundary conditions are imposed. Then let A = Ai U A 2 , on which our Hamiltonian will be 
defined. 

We next introduce fermion operators. As usual, we denote by c x>(J and c x a the annihi- 
lation and creation operators, respectively, of an electron with spin a =|, [ at site x G A. 
They satisfy the canonical anticommutation relations, 

{<V, c y , T } = {cl a , 4 iT } = 0, (2.1) 

and 

{ C i,cr) C y,r} = 3x,y3a,T (2-2) 

for any x, y G A and a, r =f, J,. We denote by N e the electron number and by $o a state with- 
out electrons. An A^-electron state can be constructed by operating N e creation operators 

We also introduce other fermion operators which play an essential role in our model. For 
each x G A, let y> x = ((p x {y)) y <=A be a vector whose components are given by [16] 

, N f 1 if \y — x\ = 1/2, /n oN 

^) = { otherwise ' ^ 

(see Fig. [2]) and define 

= E <Px{y)Cy,a = C x _l ja + C x+ i a . (2.4) 
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Figure 2: The components of the vectors (f x and cp x . 

This operator corresponds to a single-electron state localized at a bond of nearest neighbour 
sites in a chain. Since the set of all vectors ip x is not orthonormal, the 6-operators do not 
satisfy the canonical anticommutation relations (in fact, it is easy to see that {b x CT , b y ^ a } = 1 
if \x — y\ = 1), so that we introduce dual operators b x>a as follows. First recall that we adopt 
the periodic boundary conditions, and thus, for each x G Ai, sites y G A 2 can be written as 
y = x + 1/2 + n(x, y) where n(x, y) is an integer in < n(x, y) < L — 1. Then, for x G Aj, 
let ^ = (ip x (y))y £ \ be a vector whose components are given by 

!l/2 if y G A 2 and n(x, is even 
— 1/2 if y G A 2 and n(x, y) is odd (2.5) 
ifyGA 1; 

(see Fig. [2]) and, for x G A 2 , let (£ x be a vector obtained by (p x (y) = fix-i/iiv ~ 1/2)- It is 
easy to see that the inner product of £p x and <p y satisfies 

(<Px> Vy) = J2 <Pt( Z )<Pv( Z ) = <W ( 2 - 6 ) 

ZGA 

By using <p x , we define 

b x ,cr = Vx(y)Cy,a- (2.7) 
y&A 

The 6-operators are dual to the 6-operators in the sense 

{*3U> = (£x> ¥> tf )<W = 4 ; A,r. (2.8) 

The anticommutation relation (12.81) implies that the set {b x ^ojxeA is linearly indepen- 
dent and so is the set {b x ^olseA- Furthermore, (12. 8p implies that the c-operators can be 
represented in terms of the 6-operators as 

j/GA t/£A 
2/eA 
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Let us define the number operators n X}(7 with a =f, j and the spin operators S x>a with 
a = 1,2,3 for the c-operators by 

Tlx,(j c xa c xcn (2.10) 

Sx,i = - (4,T C *4 + 4 a c XiT ) , (2.11) 



S x ,3 



2 

2j (4,t Cx >.l 4,j. c £,t) ' ( 2 - 12 ) 
2(4,^-4,1^4)- ( 2 - 13 ) 



We also define 



n x = n xA + n xd . (2.14) 

The number operators n x a and n b x , and the spin operators S x a for the 6-operators are defined 
similarly by using b x>a in place of c x>a in (12.101) - f)2.14p . 

By using the operators introduced as above, we define the Hamiltonian H as follows: 

#hop = «EE bt >a b x , a , (2.15) 
xeA o-=n 

H^,J = -JY. A (^T + S x -S x y (2.16) 

H iri t,u = U^2n X) in x ,i, (2-17) 

H = H hop + H intjJ + H inttU , (2-18) 
where t, J > and £/ > 0. By using the c-operators, -f^op is rewritten as 

ffho P = t E £ 4,^ + 2* 5>*. (2.19) 

z.j/GA ct=T,| xGA 
|a-y|=l 

which represents the motion of electrons in the two chains. The Hamiltonian H int jj represents 
the repulsive interaction between two electrons with the opposite spins at the same site 
and i/hop + H int ^u reduces to the usual Hubbard Hamiltonian. In our model, owing to 
the Hamiltonian H- mt) j, an electron at a bond in one chain and an electron at a site in 
the other chain feel the attractive and ferromagnetic interaction. We stress that all the 
interaction terms considered here are of short-range. It is also noted that the Hamiltonian 
H conserves the electron number and possesses spin-rotation symmetry. The occurrence of 
ferromagnetism and electron-pair condensation in our model is thus not trivial at all. 
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Let [G X! y\ X! y e \ be an antisymmetric matrix whose elements are given by 

G x ,y 



' -1/2 if x G Ai, y G A 2 and \x - y\ — 1/2 

1/2 if x G A 2 , y G Ai and \x - y\ = 1/2 (2.20) 



otherwise, 

and define ( aT and ( UT , corresponding to pair states with spin-1, by 

Car = E G Xty b x ,J>y, T (2.21) 
x,yeA 

and 

Cut = E/ Gx,y C x,a c y,T-> (2.22) 
x,y£A 

respectively. 

Our main results are summarized as follows: 

Proposition 2.1 Assume that the electron number N e is an even integer in 2 < N e < 
\A\, and consider the Hamiltonian H with 2t — J and U > 0. Then the fullypolarized 
ferromagnetic pairing state, 

N c 

$g = (C{ T ) 2 $o, (2.23) 

is the ground state of H . Furthermore the ground state is unique up to degeneracy due to 
spin-rotation symmetry. 

Proposition 2.2 Consider the Hamiltonian H with 2t = J and U > 0, and take a sequence 
of the ground states $g of H for even N c and A such that the electron filling factor N c / (2|A|) 
converges to v in < v < 1/2. Let A — Cii/L and g k = 2cos(fc/2). Then we have 

M^) = , lim { % G) > -= (1 - 2v) I{y) (2.24) 

with 

Ily) = 2 (1- ( X [\g k \ 2 < e(v)]\g k \ 2 dk) , (2.25) 

\2TT J\k\<TT J 

where x[X] equals 1 if X is true and zero otherwise, and e{v) is determined by 

^K^i^ x[k,2 - e( "H' (2 - 26) 
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3 Proof 

Proof of Proposition 12.11 In the proof we fix the electron number iV e to an even integer 
in 2 < iV c < |A|. 

By some straightforward calculations, we can rewrite H int J as 

z xeA 

Noting the assumption It = J and the anticommutation relations (12.11) . (12. 2p and {b x a , c x ^ a } = 
0, we obtain 

H = H intyJ + H mt ^u (3.2) 

with 

#int,J = it I](4,T C ^T + 4,| c ^l)(4,T^t + ( 3 - 3 ) 
Since Hi ntt j and H irit y are sums of positive semi-definite operators 

2 (4,t c *>T + 4a c *,i) (4,T 6 *.T + 4,i 6 *.i) (3- 4 ) 

and 

Un xA n xA = Uc\^c xA c xA c xA , (3.5) 

respectively, the energy eigenvalues of #i n t,j + H- int u are greater than or equal to zero. 
Therefore, one can conclude that a zero-energy state $ of both of these Hamiltonians, which 
satisfies 

(4Ar + 4Ai) $ = ° (3-6) 

and 

Cx,iCz,T$ = (3.7) 

for all x G A (if it exists), is a ground state. We shall show that the state $g is indeed a 
zero-energy state of both H int J and H- int)Ul and thus is a ground state of H. 

Since there is no electron with down-spin in <3>g, the operation of c x ^b x ^ + c x ^b x> i on $g 
reduces to c^fe^^G- Here, by using (12. 8ft and (12. 9ft . we have 

4,t^.tC|t = 4,T {( — 1)' +1 4,T + C||^,t} = C||4,T^>T (3-8) 

for x e Ai with 1 = 1,2. Therefore, we have c x ^b x ^&Q = 0, from which ffi ntj j$G = follows. 
Furthermore, noting again that there is no electron with down-spin in $g, one immediately 
finds H int jjQ G = 0. We conclude that $g is a zero-energy state of H int J + H int>u . 
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In what follows, we shall show that any zero-energy state of H- m t,j + Hmt,u must be 
expanded in terms of $g an d its SU{2) rotations. 

Assume that $ is a zero-energy state of H and thus satisfies (13.61) and (13.71) for all x G A. 
We firstly show that <3? can be written as p2] 



*= e E/(^^T)fnOfn4j^ 

CcAi <T,T Wc* / Wc / 



(3.9) 



|C|=AT e /2 



with coefficients /, where cr and r are short hands of spin configurations (cr x )xec and (r x ) x6 c r 
with a x , t x =|, I, respectively, and the sum J2(t,t runs over an possible spin configurations. 

To prove the above claim, we begin with preparing basis states for iV e -electron states. 
Noting the linear independence of 6-operators and the relation (12.91) . we find that all states 
of the form 




* a (s T ,B i; c7 t) c7 4 ) = n Ka n b ii n 4, T n <l ^ ( 3 - iq ) 



with subsets B a , C a of A 2 such that Yjo=i \(\B a \ + \C a \) = -W e form complete basis states for 
the A^-electron Hilbert space. By using these basis states we expand $ as 

* = E E f2(B h B l ;C h C l )<S> 2 (B h B l ;C h C l ), (3.11) 

B T ,B|CA 2 C* T ,C* i cA 2 

with coefficients f 2 , where the sum runs over all possible subsets of A 2 such that J2a='\,i(\B ( j\ + 
| Co- 1) = N e . (For notational simplicity, here and in the rest of the proof we do not explicitly 
write this restriction in the expressions.) 

Let us derive the condition on f 2 which follows from (13.61) and (13.71) . We choose a site 
z e A 2 and B' a , C a C A 2 with a =f, | satisfying z £ B' a and E^fl-B^l + \C a \) = N e -2. Then, 
from ( 13. 6p . we have that 



n c ^i) ( n c -,t ) k-^zA^ z m + = o, (3.i2) 




which leads to f 2 {B\ U {z},B{ U {«}; C T , Q) = for z <£ C a . Therefore, if z £ C T n C h 
then f 2 {B\ U {z},B[ U {z};C h C±) = or equivalently / 2 (£ T , C T , = for 5 T and 
B± such that 2 G B^ D f?j,. On the other hand, by examining (13.71) for sites in A 2 , we also 
obtain f 2 (B^, B±; Cj, CjJ = for C-j- and Cj_ with C|nC| 7^ 0. As a result we find that 
f 2 (B h B^ C h Cjj = unless 5 T n 5 X = as well as C T n Cj_ = is satisfied. 
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We furthermore examine (13.61) . Fix Cj and C± with C| U C± 7^ A2 and choose a site 
z G A 2 which is not included in Cf U Cj. Using the condition (13.61) for z, we have that 



n ^,1) ( n ft.,t) f n ^) [ n ^ <v(4at + 4,^)$ = o (3.13) 

eB[ J \xeE'^ J \x€C L J \x£C T / 



1. From this we obtain 



for arbitrary B\, B[ C A 2 such that Ea(\K\ + \ C A) = N c 
f 2 {B\ U {z}, B{; C h C x ) = and f 2 {B\, B{ U {z};C t ,C;) = for 2 ^ C T U C x , i.e., that 
/ 2 (S|, Sjj C|, C|) = if there exists a site z such that z £ B| U B| and z ^ Cf U C^. 
Therefore, we conclude that ^(-Bf ^T> ^i) = unless all the conditions Cf D C± = 0, 
5 T n B i = and (5 T U B|) C (C T U Q) are satisfied. 

Taking into account the above result and noting ( 12.7ft and (12.91) . we expand $ as 



E Ex 

B T ,Bj,CAi C T ,Cj,CAi 



EW>EI^I 



h{B h B i; C h CJ^Bi, B i; C h Q), (3.14) 



with coefficients fi, where we define $1 as in ( 13.101) . replacing subsets C a and B a of 
A 2 with those of Ai. Then, repeating the same argument as above, we conclude that 
fx{B h B^C h C{) = unless C T D Q = 0, £? T n B t = and (5 T U B\) = (C T U Q). 
We thus reach the desired expression (13.91) of a zero-energy state $. 
We shall next show that the coefficients / in (13.91) satisfy 



/(0,t)=/(C";<t',t') 

i)- 

For z G A 2 and C C Ai a straightforward calculation yields 



(3.15) 



'n*u)(n4j*o 



E ^/(^W^)^ 11 ^^^] (xWv =T>°V =1} - xWv =i^v' =TD 



x 



n &u)(n4.r.W 

^GC\{j/,y'} / Wc / 



(3.16) 



Here sgn[- ■ •] is a sign factor arising from exchanges of the fermion operators. By definition, 
for any y, y' in C C Ai, (p y {z)(p y i{z) with z G A 2 is non-vanishing. We thus have from (13.71) 
and (13 . 161) that /(C;<x, t) = f(C;cr[y <-> J/'],t), where <x[y <-> y'] is a spin configuration 
obtained by switching a y with cy. 
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It is also easy to see that 

(4Ar + 4Ai)(n&u) fniJ 

\xeC / \xeC I 
= sgn[z; C] (xWz =T, r g =1} - x[<r z =1, r z =T]) 

xiAl ( II &U ) f II 4,J (3.17) 

\xec\{z} } \xec\{z} ) 

for z G C C Ai. We thus have from (J3l| and fl3~TTj) that /(C; cr, r) = /(C; cr[ 2 ], r[z]) where 
<t[z] and t[z] are obtained by switching a z and t z in cr and r. 

So far we have shown (13.151) with C = C . We shall complete the proof of the claim by 
examining (13. 6j) with z G A 2 . 

By using (12 .4p and (12.91) . (13.61) is rewritten as 

£ ^KG/X^',! + ^ 4 ^'4) $ = °- ( 3 - 18 ) 

Let ?/,?/ be a pair of sites with \y — y'\ — 1 in Ai. For this pair, there exists a site z in 
A 2 such that ip z (y)ip z (y f ) ^ 0. Let D be an arbitrary subset of Ai such that y,y' ^ D and 
|-D| = A e /2 — 1. Then, noting (I3.18p . we obtain from 



n l 

\xeD 



J) ( II ^)«T 6 *.T + c U 6 -.l) $ ( 3 - 19 ) 

/ \xeDU{y,y'} ) 

o f n f n ^ ] + ^,v^-v) $ = ° ( 3 - 2 °) 

VrcGD / \xe-DU{y,y'} / 



= (Pz(y)<Pz(V 

that 

/(£> U {y}; cr, r) = /(£> U {j/'}; cr', r') (3.21) 

where a y = t y >, r y = t y , a y > = i y , r y i = i y i and a x = a' x = v x , r x = r' x = i x for x G D. Let C 
be an arbitrary subset of Ai with \C\ — N e /2 and let C y ^ y i be a subset which is obtained by 
removing y from and adding y' to C (we assume y G C and y' ^ C). Then the above result 
implies that (13. 15[) holds for C" = C y ^ y i if y and j/' satisfy \y — y'\ — 1. It is easy to check 
that, for arbitrary subsets C and C" of A x with |C| = |C"|, there exist pairs ({?//, 2/;})j=i with 
| j/, - y[\ = 1 such that C 1 = C yi ^ y[ ,C 2 = C x m ^---,C = C^ yln . This completes the 
proof of the desired relation (13.151) . 

We write f(M) for /(C; <r, t) with + T x) = M. By using fl3TT5|) . (ET5j) is 

rewritten as 

* = £ E E/Wx[E(^ + r.) = M](nOfn4j^o 

CcAi M=-N e /2 <T,T xeAi VieC / VreC / 



./ 

|C|=JV e /2 
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N *l 2 I^-AT ( \ ( \ 

= e e /'(M)(^ t ) 2 ns,t n< t *o 

CcAi M=~N c /2 VreC I \x&C / 

\C\=N e /2 

= E f"(M)(S- ot ) 2 $ G (3.22) 

M=-N c /2 

where 

f(M) = W/2 + M ^ /2 - M " /(M). (3.23) 

/"(Af) = (-!)*(*-') _J_/'(Af), (3.24) 

and the total spin lowering operator S^ ot is defined as S^t = ExeA c l,j. c x,T- To get the final 
expression we used 

Ctt = E Gx,ybx,tKj + E G x ,yb x ^b y ^ = E c x,T^,T- (3.25) 

xeAi,2/£A2 xeA 2 ,j/eAi iGAi 

This completes the proof of Proposition 12. 11 

Proof of Proposition I2T21 As in the case of Ctt; we have 

Ctt = E ^,T C *,T- (3.26) 
xeAi 



The momentum representation of is 



where 



Ctt = E fl , feC(2,fc),TC(i,-fe),t, (3.27) 
fce/c 



&(l,k),a = ^=E e " ifcXc ^' ( 3 - 28 ) 



(3.29) 



with / = 1, 2 and 



/v' [k = —n 



L-l L-l 



<n< — — \ . (3.30) 
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By using (13.251) and (13.261) . we have 

*g oc e f n biA,i) $ o oc e f n %,A,r) ( n ) ^ 

CcAi \xgC I CcAi ViGC / \xgAi / 

\C\=N c /2 \C\=N c /2 

<* E ( II &*, T c,, T )$FOc(Cn)^$F (3.31) 

CcAi \xeAi\C / 

|C|=JV e /2 

where $f = (ilxeA c l,t) and iVh = |A| — jV e , which is the number of holes counted from 
half-filling of electrons. We use the above representation of $g with in (13.271) to estimate 
a lower bound on 

It is easy to check the anticommutation relation {ct ^ a , C(jr ; k'),a'} = 8i t i>5k tk '8a,a' , and by 
using this relation we have 

CttCtt = E M 2 -EE \9k\ 2 c { i, k )A-c{i, k) ^ + CttCV ( 3 - 32 ) 

fcG/C fcG/C Z=l,2 

We furthermore have 

C(i,fc),Tc[i ;fc ) >T Ctt = ^-fcC( 2 ,-fc),tC(i,fc),t + CnC(i,fc),tc[i,fc),T' ( 3 - 33 ) 

which leads to 

£ (i,fc).T £ (i,*),T ^Tt)" 1 $ f = mg-_ fe c (2 _ fe)iT c (1)fe))T (Ctt)" 1 " 1 $ f (3.34) 
for integers m and thus 

(C n )-" m 6(i.*),TcJi,*),T ^t)™ $ f = ^-c (1 , fc)iT cf lifc)iT (C TT )~ *f. (3-35) 
The same relation holds for C(2,fe),fcL ^ j. We thus obtain 
($ G ,AAt$ G ) _ N±J_^, |2 



1 fN h \ 1 ($ G , C(i,fc),TC(/,fc),T $G > 

2 It - X J £ S i»i 



iV h 1 „ „ |2 ( $ G,c{w,T a (W 5 T $ G) 

> ~^J2 2. 1^1 ($ $-\ 

^ n keKl=l,2 \^G)^G/ 
I A I — N 2 ^ 

> W. (3-36) 

where we arranged the elements in K as k(l), k(2), . . . , k(L) so that \gk( m )\ 2 < |<7fc(m')| 2 if 
m < m'. Taking the limit |A|, N e ^ oo with iV e /(2]A|) converging to u completes the proof. 
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4 Remarks 



In this section we remark some aspects of our model and its ground states. 

The state (S^ oit )~2'~ M $G is the representative of the ground state in the subspace where 
the third component of the total spin is M. If we use the expression ( 13.311) of $g> this 

becomes (Ctt)^^m where $ M = (5' t ~ t )^ _M $F- This implies that the ground state can be 
regarded as a hole-condensation state in which all holes form the spin-triplet pair in the back 
ground of the fullypolarized ferromagnetic state. 

By using the so-called d vector d = (di, e^, d^), one can obtain a useful representation of 
a pair operator [18]. In our case, we define 

ci= E E F Zid)liXr (4- 1 ) 

x,y£Ai <t,t=U 



with 




x,y(d) 




(4.2) 



The pair operator £L corresponds to the case of the complex d vector d = (—1/2, — i/2,0), 
and thus the ground state $g is a nonunitary spin-triplet pairing state similar to the one 
describing the Ax-phase of superfluid 3 He in the magnetic field [18]. It is noted that, unlike 
the case of Ax-phase realized in the magnetic field, the states 

t ^ 2 



Q) 2 $ o (4.3) 

with any d obtained by rotating (—1/2, — i/2, 0) are the ground states of H, since our Hamil- 
tonian H has spin-rotation symmetry. The direction of the magnetic moment is given by 
the cross product id x d* where * means the complex conjugation. 

Let B be the collection of pairs {x, y} of sites such that x G A 1; y 6 A 2 and \x — y\ = 1/2. 
Then, we have from (13.261) that 



{x,y}ee 



Therefore, the order parameter fi(i>) measures the long-range correlation between local spin- 
1 electron pairs, and Proposition 12.21 implies that there exists this long-range correlation in 
the ground state of H in the thermodynamic limit. 

With respect to the ground state $g expectation values of operators constituted of either 
annihilation or creation operators, such as A and Cf2,k),iCn-k),i, are always zero since there 
are exactly N e electrons in $g- m order to obtain a particle number symmetry breaking 
ground state, which is usually discussed in mean field approximations, we need to form a 
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linear combination of $g with different electron numbers. Here let us consider the Bardeen- 
Cooper-Schrieffer type state 



$ G = II I 1 + 9kC{2,k),-\kh-m) $F- (4.5) 

fce/c 

It is noted that the projection of $ G onto the Hilbert space with the fixed electron number 
is proportional to $ G an d thus $' G attains the ground state energy of H . Then, for $ G the 
expectation values of pair annihilation and creation operators are calculated as 

($G» C(2,fc),tC(l,- fc ), T $ G ) = (go, jl,-k)Am^G) 9k ( fi x 

($'g, *g> (&q,#q) 1 + s2' K ' J 

which is finite for — tt < k < tt. 

When we set U = 0, the model with 2t = J has degenerate ground states and does not 
exhibit ferro magnetism. In fact, states of the form 

(cy""^)"" (c| t ) s " *. (4.7) 

with non-negative integers N^, and are zero-energy states of H mt j, since we have 

(4,T 6 *>T + 4Ai)Sr = Cir(4,t 6 *,T + 4,1^4) ( 4 -8) 

for any x G A and a, r =\ , j, as in (13.81) . The on-site repulsion removes this degeneracy and 
generates the unique ferromagnetic pairing ground state. 

We can extent the present model, which has spin-rotation symmetry, to anisotropic spin- 
interaction cases. Consider the Hamiltonian 

H a >P = H hop + Hggj + H intiU (4.9) 

where H^f j is given by 



H-Z = -^E|^ + PS b xA S xA + (3S b x>2 S Xj2 + (a-aP + (3)S b X;3 S x , 3 j (4.10) 

with a = ±1 and /3 in < /3 < 1. It is noted that H ±1,1 becomes H. 

The following results are obtained for the Hamiltonian H a, @. In the case of a = 1, < 
(3 < 1, 2t = J, U > 0, for even iV e in 2 < ]V e < |A|, there exist exactly two ground states of 
H 01 ' 13 , which are given by 

N c 

*g = (Cl) 2 $0 (4.11) 



b 
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with a =|, j. In this case, therefore, the ground states of H a ^ exhibit ferromagnetism and 
condensation of the spin-1 electron pairs whose spins point in the same direction parallel to 
the third component axis. In the case of a — — 1, < /3 < 1, 2t = J, U = 0, the ground 
state of H a, @ is unique for even N e in 2 < N e < 2|A|, and is given by 

$ G = (C||) 2 $o. (4.12) 

In this case the ground state exhibits condensation of the spin-1 electron pairs whose spins 
point in the plane perpendicular to the third component axis, but does not exhibit ferro- 
magnetism. One can prove these results in the same way as in section [3j by noting the 
relation 

= pH int>J + (1 - (3)H« tJ + H^u (4.13) 

with 

^int,J = IT b X(T c Xja(T c xa(T b XtCr . (4-14) 

1 xeA<r=U 

In this paper, we treat only the case of 2t — J where the exact spin-1 pairing ground 
states can be constructed. Apart from 2t = J we have the following exact results. When we 
set J = 0, the model becomes the decoupled two Hubbard chains, and thus the ground state 
of H is non-magnetic. In another limit of U —>■ oo, where double occupancies at sites are 
forbidden, it can be proved, by using the Perron- Frobenius theorem [19], that the ground 
state of H for 2t < J, < N c < |A| with open boundary conditions is ferromagnetic. So we 
expect that the ground state phase diagram of our model has a rich structure depending on 
J, U, and the electron filling factor v (and also the anisotropy parameters a and j3 in the 
anisotropic cases). It is desirable to describe the phase diagram in detail, but it is beyond 
of our scope at the present time. We leave this problem for future study. 



5 Higher Dimensional Cases 

Let L be a positive odd integer and define K\ = [0, L — l] d D Z d . Let a be a vector whose 
all components are 1/2, and define A 2 = {x + a | x £ Ax} and A = Ai U A 2 . We impose 
periodic boundary conditions in all directions on A. For each x e A, let cp x = (<p x (y)) y <=A 
be a vector whose components are given by 



<p*(v) = \l if 1^*1 = 1-1 = ^ ({U) 

|^ otherwise. 

Recalling that the periodic boundary conditions are adopted, one notices that for each 
x E Ai, sites y E A 2 can be written as y = a + x + Yl'i=i n i{ x ^y)^i where Si is the unit 
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vector along the /-axis, and ni(x, y) with I = 1, . . . ,d are integers in < ni(x, y) < L — 1. 
Then, for x G Ax, let = ((p x (y)) y€ \ be a vector whose components are given by 

!l/2 d if y G A 2 and Yjf =1 "■«(#, 2/) is even 
-l/2 d if 2/ 6 A 2 and £f =1 ni(x, y) is odd (5.2) 
if y e Ax , 

and for cc G A 2 , let y?^ be a vector obtained by <fi x (y) = <p x -a(y — a )- Let [G^y^yeA be an 
antisymmetric matrix whose elements are given by 

!— 1/2 if x G A 1; y G A 2 and |aj — j/| = |a| 
1/2 if x G A 2 , y G Ax and jcc - y\ = \a\ (5.3) 
otherwise. 

By using ip x , cp x and [G^.y^ygA introduced as above, we define b- and 6-operators and 
pair operators C CTT and £ aT , and then define the Hamiltonian H, as in section [2j For this 
(i-dimensional Hamiltonian we can obtain the same results as in Proposition 12.11 and 12.21 
where I{y) is replaced with 

I i (y) = 2( 7 ^ s f ■■■[ x[|^| 2 <6(z/)]|^ fc | 2 dA;i---d^y (5.4) 

Here g^ = 2 d Y\f =1 cos(fc;/2) and e(y) is determined by 

"=UtTu[ X[\g k \ 2 <<iy)}dk 1 ---dk d ). (5.5) 

2 \(2n) d J\ki\<n ^|/fc d |<7r / 

We can also obtain the <i-dimensional anisotropic spin-interaction model as in section HI 
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